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g&bstract 

o : 

^lt is shown that the existence of hot accretion spots at the surface of T Tau prac- 
tically has no effect on the accuracy of estimation of its magnetic field strength at 
^■photospheric level. 

t-h We also found that one can interpret results of T Tau's photospheric magnetic 
^4ield strength measurements carried out via different methods in the frame of the 
following alternative: 1) if T Tau's inclination angle i < 10°, then magnetic field 
^bf the star may be dipolar with the angle between rotational and magnetic axes is 
oiiear 85°; 2) if it will be found (e.g. from interferometric observations) that i > 10°, 

i/^hen magnetic field of T Tau is essentially non-dipolar or/and non-stationary. 

O ■ 

^Introduction 

o ■ 

Tauri stars are low mass (M < 2M ), young (t < 10 Myr) pre-main sequence 
j^tars. T Tauri itself belongs to subclass of classical TTS, what means that its activity 
j^s the result of (protoplanetary) disc accretion onto central young star. The problem 
^J)f strength and geometry of TTS's magnetic field is one of the most important in 
the physics of young stars because it is widely belived that magnetic field determines 
the character of TTS's activity and controls evolution of stellar angular momentum. 
At the moment the strength of magnetic field is known for less than ten TTS due 
to difficulties of magnetic field measurements in these objects. 

Basing on the fact that magnetic field changes equivalent widths of spectral lines 
depending on the value of their Lande factor and line splitting pattern, Guenther 
et al. (1999) and Johns-Krull et al. (1999) found that avereged surface magnetic 
field strength of T Tau B > 2 kG. On the other hand, Smirnov et al. (2003, 2004) 
found from specrtopolarimetric observations that averaged longitudinal (i.e. line of 
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site oriented) component of T Tau's surface magnetic field B\\ was ~ +150 ±50 G in 
1996 and 2000 yr, and ~ +50 ± 30 G in 2003 yr, i.e. more than order of magnitude 
less than 5-value. Possible reasons of so strong discrepancy is considered below. 

The influence of accretion onto accuracy of estimation of 
observed B and B\\ values 

Calculations of Romanova et al. (2003) demonstrated that when disk accretion 
occurs onto a star with dipolar magnetic field gas predominantly falls near stellar 
magnetic poles, i.e. onto regions where magnetic field strength is maximal. Accreted 
matter heats stellar surface, so if contribution of "hot" spots into formation of some 
spectral line is negligible, then average B and B\\ values, derived from observations 
of this line, will be less than in the case of star without hot spots. Let estimate this 
effect in a quantitave way. 

Effective temperature T e f of the star in the accretion zone can be estimated from 
the energy conservation equation: 

aT? + ^-~ SI*, (1) 

where T* is an effective temperature of the star outside the accretion zone; p and 
V are density and infall velocity of accreted matter respectively. Spectra of regions 
with T > T* differ from spectrum of undisturbed photosphere. To take this effect 
into account we will suppose in the following that spectral lines of interest do not 
originate in these regions at all. 

We adopted p = p(9, <p) and V = V(9, p) distributions of infall gas parameters 
along stellar surface - 9 and p are latitude and longitude respectively - from Ro- 
manova et al. (2003) calculations and used them to estimate T e f{9, tp) distribution 
from Eq.(l), assuming = 5250 K (White & Ghez, 2001). Three dimensional 
magnetohydrodynamic calculations of Romanova et al. (2003) were carried out for 
different values of an angle a between magnetic and rotational axes. 1 Temperature 
distribution shown in Fig.l was calculated with a = 75°. More precisely, lines of 
AT (9, p) = T e f — T* = const, that characterize excess heating of stellar surface in 
the accretion zone are shown in this figure. Similar distributions for other values 
of a-parameter can be constructed from data presented in Romanova et al. (2004) 
paper. 

1 In addition, Romanova et al. assumed that polar field strcnght of the star Bo is ~ 2.1 kG and accretion rate 
M ac ~ 2 • f(T 7 M Q /yr. 
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B\\ and B values that one derives from observation of photospheric lines results 
from averaging of respective local values over visible hemisphere of the star with 
radius R*: 

B = ^2 hB(e,tp)coBpdS, B W = ^2 f Sm B(e,tp) cosficoB-ydS. (2) 

Consider spherical system of coordinate with the origin in the center of the star, 
9 = axis of which coinsides with stellar rotational axis, characterized by unity u 
vector. Let <p = plane of this system be the plane that passes throught vector u 
and unity vector 1 oriented to the Earth. Let also n be a unity vector normal to 
stellar surface at given point and B be a vector of magnetic field at the same point. 
Then in Eq.(2) R~ 2 dS = sin 6 d6 dip; (3 is the angle between n and 1 vectors; 7 is 
the angle between 1 and B, such as: 

B 

cos P = sin i sin 6 cos p + cos % cos 6, cos 7 = 1—, 

B 

where i is the angle between stellar rotational axis and the line of site, i.e. between 
to and 1 vectors. 

B(#, ip) vector will vary with time t due to stellar rotation. Let assume that 
magnetic field of the star is dipolar, i.e. 

B 3n(/7n) — /2 
B~o = 2 ' 

where Bq is polar magnetic field strength and /2 is the unity vector oriented along 
stellar magnetic axis. Then to discribe B(t) dependence in each point of stellar 
surface one needs only one parameter - phase ip (0 < ip < 1), that characterise rota- 
tion of stellar magnetic pole relative to <p = plane. Accretion flow disturbs initial 
(dipolar) magnetic field of the star, but it follows from calculations of Romanova et 
al., (2003) that the difference between initial and disturbed fields near stellar surface 
does not exceed 1 %. 

To estimate the effect of hot spots existence on the average values of total mag- 
netic field strength and its longitudinal component we calculated "observed" values 
B\\ and B by means of Eq.(2) for three cases: 1) neglecting with existence of hot 
spots, i.e. calculating integrals over all visible hemisphere of the star; 2) excluding 
regions of stellar surface (accretion zone) where AT > 200 K; 3) excluding regions 
where AT > 1000 K; 

Variations of B^/Bq value with phase of stellar rotational period ip for these three 
cases are shown in Fig. 2 as an example of our calculations carried out with a = 75° 
and four values of inclination angle % : 10°, 30°, 60° and 90°. Analogous variations 
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B(ip)/Bo for the same values of a and i parameters are shown in Fig. 3. It can be 
seen from these figures that the difference between curves does not exceed 10 — 15 % 
at any value of i. Similar result was found at other values of a parameter. It means 
that errors in measured values of B and B\\ due to neglecting of existence of hot 
spots is less than 15 %, what is less than modern accuracy of measurements of these 
values in TTSs. It is the consequence of relatively small (< 10%) area of stellar 
surface occupied by accretion zone, where the main portion of accretion energy is 
liberated (Romanova et al. 2003; 2004). 

Estimation of parameters of T Tau's magnetic field under 
assumption that it is dipolar 

Let assume that global magnetic field of T Tau is dipolar and stationary. Then 
more than order of magnitude difference between observed values of B\\ and B is 
due to proximity of angle 5 between magnetic axis and line of sight to 90° at the 
moment of observations. (At 5 = 90° longitudinal field components that parallel 
and antiparallel to the line of sighte are mutially compensated and average value of 
Bu should be exactly equal to zero.) Consider this question in a quantitative way 
basing on the results of existing measurement. 

To discribe dipolar magnetic field one needs three parameters: Bq, i and a. For 
given set of these parameters observed (average) values of B and Bu will vary with 
phase ip, e.g. as shown in Fig. 2 and Fig. 3. Let suppose that probability p to detect 
value B in the range B obs ± 3 a at some moment is equal to the lenght of phase 
interval, inside which | B — B obs |< 3d, where a is the error of measurements. 
Probability to detect value B\\ in the range | B\\ — B^ bs |< 3d is determined in a 
similar way. 

Guenther et al. (1999) and Johns-Krull et al. (1999) found B = 2.3 ± 0.15 kG 
and the following values of Bn was found in our observations: +150 ± 50 in 1996 
yr and 2002 yr as well as +50 ± 30 G in 2003 yr. If Bq and a parameters did not 
vary with time (magnetic field is stationary) then the difference in observed values 
of B and B\\ is the consequence of observations at different phases ip. Then one can 
find for any set of three parameters Bq, i and a probability P = Pi ■ P2 ■ Ps ■ Pa ■ Pb 
to detect two 5-values and three Bu that was indeed found from observations. 

Such calculations were carried out in the following range of parameters: Bq from 
1.5 to 5 kG, 0° < a < 90° and 0° < i < 60°. We took into account that % = 29 ±l% 
according to Akeson et al. (2002) when adopting the range for i. Calculations 
indicated that: 
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1) Probability P differs from zero only if 3kG< Bq < 4.3 kG at any values 
of a and % from the range of interest. This result looks natural if to take into 
account that surface magnetic field varies between Bq (at magnetic pole) and Bq/2 
(at magnetic equator) in the case of dipol, so average B-value should be inside the 
range (B /2,B ). 

2) P = P(a, i) function is practically independent on values of Bq from allowable 
range. As an axample we plot in Fig. 4 P = P(a, i) function calculated at Bq = 3.5 
kG. (Lines of equal probability consist of segmnets of stright lines because P(a, i) 
was calculated at discrete values of a and i angles with 5° step.) 

As follows from Fig. 4 at a < 80° and/or % > 10° probability to find observed 
values of B and B\\ is less than 0.2, so one can consider this possibility as a low 
probability one. (Note by the way that this conclusion is practically independent 
on the real value of T Tau's Bq.) Two conclusions follow from this statement: 

1) If future observations (e.g. with VLTI interferometer) will demonstrate that 
inclination angle of T Tau is larger than 10°, then existing results of B and B\\ 
measurements mean that magnetic field of the star is essentially non-dipolar near 
its surface. Note that there are some theoretical reasons to belive that configuration 
of magnetic fields of TTSs can be close to quadrupol (Dudorov, 1995). 

2) If it will be found that % < 10°, then existing measurements of magnetic field 
in T Tau do not contradict to the hypotesis that it has dipolar magnet field with 
dipol axis inclined to the rotational axis at a ~ 85°. 

Note once more that these statements are valid only if magnetic field of T Tau 
is stationary, otherwise there is no sence to compare values of B and B\\ that was 
observed non-simultaneously. 

Conclusion 

Calculations of Romanova et al. (2003) indicate that in the case of disc accretion 
onto inclined magnetic dipol, with parameters Bq and M ac typical to classical TTS, 
the difference between resulting and undisturbed magnetic field near stellar surface 
is negligible. On the other hand we found that in the case of dipolar magnetic field 
the existence of hot accretion spots at stellar surface practically has no effect on the 
accuracy of measurements of average B and B\\ values. 

It means that comparison of observed B and B\\ values at photospheric level 
could give information to understand if TTS's magnetic field is dipolar or has more 
complicate geometry. We demonstrated above that in the case of T Tau one can 
agree observed values of B and B\\ in the frame of the following alternative: magnetic 
field of the star is eccentially non-dipolar or it is dipolar, but dipol axis is inclined 
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to stellar rotational axis at ~ 85°. It will be possible to discriminate between these 
possibilities when future (interferometric) observations will answer the question: 
does inclination i of stellar rotation axis to the line of sight exceeds 10° or not. 

It would be more precisely to say that such alternative follows from observations 
only if magnetic field of T Tau is stationary. If it does indeed - this question is 
interesting by itself but to answer it more observations are necessary. 

This work was supported by grants RFBR 02-02-16070, INTAS 03-51-6311, NASA 
NAG5-13220, NAG5-13060, NSF AST-0307817. 
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Figure 1: Distribution of excess temperature AT on the surface of a star that has dipolar magnetic 
field with the angle between magnetic and rotational axes a = 75° (see text for details). The 
innermost countour respects to AT = 1000 K. Longitude and latitude (in degrees) are plotted along 
X and Y axes respectively. 
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Figure 2: Longitudinal component B\\ of stellar magnetic field strength as function of phase of 
rotational period ip for a = 75° and different values of inclination angle i. Dot-dashed lines respects 
to situation when all stellar surface is considered (there is no hot spots), solid lines - stellar surface 
without regions with AT > 200 K, dashed lines - stellar surface without regions with AT > 1000 K. 
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Figure 3: The same as in Fig.2 but for average magnetic field strength B. 
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Figure 4: Lines of equal probability P(a,i) = const in the case when B = 3.5 kG. See text for 
details. 
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